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We investigate jet quenching of virtual gluons and thermalization of a strongly-coupled plasma 
with a non-zero chemical potential via the gauge/gravity duality. By tracking a charged shell falling 
in an asymptotic AdSd+i background for d — 3 and d = 4, which is characterized by the AdS- 
Reissner-Nordstrom-Vaidya (AdS-RN-Vaidya) geometry, we extract a thermalization time of the 
medium with a non-zero chemical potential. In addition, we study the falling string as the holo- 
graphic dual of a virtual gluon in the AdS-RN-Vaidya spacetime. The stopping distance of the 
massless particle representing the tip of the falling string in such a spacetime could reveal the jet 
quenching of an energetic light probe traversing the medium in the presence of a chemical poten- 
tial. We find that the stopping distance decreases when the chemical potential is increased in both 
AdS-RN and AdS-RN-Vaidya spacetimes, which correspond to the thermalized and thermalizing 
media respectively. Moreover, we find that the soft gluon with an energy comparable to the ther- 
malization temperature and chemical potential in the medium travels further in the non-equilibrium 
plasma. The thermalization time obtained here by tracking a falling charged shell does not exhibit, 
generically, the same qualitative features as the one obtained studying non-local observables. This 
indicates that -holographically- the definition of thermalization time is observer dependent and 
there is no unambiguos definition. 

PACS numbers: 



I. INTRODUCTION 



The AdS/CFT correspondence, a duality between the type IIB supergravity and the Af = 4 
SU(N C ) super Yang-Mills theory (SYM) at large 't Hooft coupling in the large N c limit l|-|4|, has 
been widely studied to investigate the properties of strongly coupled systems. At finite temperature, 
the spacetime metric in the gravity dual is governed by the AdS-Schwarzschild_geometry, in which 
the temperature is characterized by the Hawking temperature of a black hole 5| • The holographic 



dual description can be regarded as an analogue of the strongly-coupled quark gluon plasma (QGP) 
generated in the relativistic heavy ion collisions. On the other hand, the QGP also carries a non- 
vanishing chemical potential. In the gravity dual, the chemical potential is encoded in the AdS- 
Reissner-Nordstrom (AdS-RN) spacetimeP, Q] which carries a no n- vanishing gauge field. 

One of the interesting aspects of the heavy ion physics is the thermalization process of the medium 
after the collisions of the two nuclei. In the gravity dual, this scenario should correspond to the grav- 



itational collapse and the formation o 



approaches in the literature. In 



a black hole. The issue was studied by implementing various 



8Hl9| the collisions of grav itational shock waves are introduced to 



mimic the colliding nuclei in the relativistic collisions. In 



20 



2j] ti 



time-dependent and boost-invariant 



metrics, associated with the plasma undergoing Bjorken expansion, were investigated and further 



generalized in 22| with a chemical potential. Moreover, the authors in 23 



24j introduced anisotropic 



and time- dependent boundary conditions, which lead to the thermalization of an anisotropic plasma. 



Alternatively, the gravitationa 



the isotropic thermalization 



25 



collapse can be characterized by a collapsing shell, which results in 



3. 



Note that, such first principle computations are generally diffi- 
cult and thus a more phenomenological approach allows us to probe the physics with more ease. 
We also note that the back grou nd metric generated by a time- varying weak dilaton field produces 



AdS-Vaidya-type background 
shell limit. In 



271 ] and provides a very good quantitative approximation in the thin- 
-434(| . various non-local operators have been studied in the AdS-Vaidya spacetime 
to probe aspects of thermalization of the medium; and generalized in the presence of a c harg ed 



32 



35 



36l | . As shown in 



37|, a 



shell where the spacetime is represented by AdS-RN-Vaidya metric 
thermalization time which is independent of the length scales of non-local operators can be extracted 
from an alternative approach. 

Another salient issue in the heavy ion collisions is the jet quenching of hard probes traversing 
the medium. Much effort has been made towards studying this phenomenon in the thermalized 
medium via holographic methods. In general, the heavy probes such as heavy quarks are assumed 



391 ] . jet quenching of heavy quarks is 



to constantly travel through the infinite medium. In [38 . 
characterized by the drag force caused by the interaction with the medium, which can be derived 
from the dynamics of a trailing string in the gravity dual. Furthermore, the jet quenching parameter, 



which encodes the momentum broadening of hard probes, can also be extracted from the computation 
of the light cone Wilson loop in the curved spacetime 



40 



42l | . However, for light pro 



3es such as light 



49(, the gravity 



quarks and gluons, they may finally dissipate in the medium. According to [43 
dual of light probes may have various candidates. The maximum stopping distance of the light 
probes can, nevertheless, be derived from the null geodesic of a massless particle falling in the dual 



geometry. There are only a hand 
the jet quenching of hard probes 



ul of studies on the influence of a non-equilibrium medium on 
Some work on the influence of thermalization on 



37 



n 3-54 



electromagentic probes can be found in 52H54|. In 37(], it is indicated that jet quenching of light 
probes with energy much greater than the temperature of the thermal medium remains unaffected 
by the non-equilibrium processes. This situation may change when the probes carry finite energy 
comparable with other soft scales of the medium. 
In this paper, we will follow the previous studies in 



35 



37| to investigate thermalization of the 



medium with a chemical potential and jet quenching of light probes in such a non-equilibrium plasma. 

Our work is organized in the following order. In section II, we analyze the AdS-RN-Vaidya metric 

in Poincare coordinate and extract a thermalization time by tracking the position of the falling shell 

in the thin-shell limit. We compare our results with the ones obtained from non-local observables 
1 

351 ] and argue that these two prescriptions can only be compared when the length scale of the 



non-local operators is roughly the size of the future horizon. In section III following 44j , we study 
jet quenching of a virtual gluon traveling in the medium, which is characterized by a double string 
falling in the gravity dual background. First we compute the stopping distance of a gluon in the 
thermalized medium with non-zero chemical potential. Then we set up the initial conditions of the 
falling string in the AdS-RN-Vaidya geometry and the matching condition when the tip of the string 
penetrates the shell. Towards the end of the section, we evaluate stopping distances for both hard 
and soft gluons in the non-equilibrium medium with a chemical potential. We conclude with a brief 
summary and discussions. Some technical details are presented in two appendices. 
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II. FALLING SHELL IN ADS-RN-VAIDYA SPACETIME 

In this section, we will extract the thermalization time of the non-equilibrium plasma with a 
non-zero chemical potential by tracking the position of the falling shell in Poincare patch of AdS- 



RN-Vaidya geometry. The approach is different 
Finkelstein (EF) coordinates as carried out in 



35 



rom studying non-local operators in Eddington- 



361 ] . Since we can define thermalization time to 



be the time when the shell almost coincides with the future horizon, it is independent of the length 
of the operators. We will observe that with increasing chemical potential, this thermalization time 
decreases. 
We begin with the AdS-RN-Vaidya 



571 ] metrics in EF coordinates with d = 3 and d — 4, respec- 



tively. For d = 3, we have 

ds 2 = — (— f(v , z)dv 2 — 2dvdz + dxf) , A v = q{v){zh — z) , (1) 

z 

where f(v, z) = 1 — m(v )z 3 + ^q(v) 2 z A . For d = 4, we have 

ds 2 = — (-f(v,z)dv 2 -2dvdz + dx 2 ) , A v = q(v)(zl - z 2 ) , (2) 
z l 

where f(v,z) = 1 — m{v)z 4: + ^q(v) 2 z 6 . Here we have set the AdS curvature radius L — 1. The 
coordinates Xj represent the d- dimensional spatial directions. Also, v denotes the EF time coordinate 
and z denotes the radial direction. The boundary is located at z = and Zh denotes the future 
horizon. In the equations above, m(v) and q(v) represent the mass and electric charge of the shell. 
The charge is related to the time component of the vector potential, which generates an R-charge 
chemical potential in the gauge theory side via 



fi = lim AJv, z) . (3) 

2->-0 



For simplicity, we are interested in the thin-shell limit of the interpolating mass function 

m(v) = — ( 1 + tanh ( — J J (4) 



2 V W 



with t»o 0. We can make the same choices as 



by taking q(v) 2 = Q 2 m(v) 4 ^ 3 and q(v)' 



Q 2 m(v) 3 / 2 for d = 3 and d = 4, respectively. The values of are determined by 

/(« > « 0j * A ) = 1 - Mzl + \^ IZ Q 2 A = , for d = 3 , (5) 
/(v > u , = 1 - Mz 4 h + 2 -M z l 2 Q 2 zl = , for d = 4 . (6) 

We assume the shell falls from the boundary at t — and thus /x is a constant since t> = t > on 
the boundary. On the other hand, the thermalization temperature is given by 



The entropy can be determined by the area of the black hole, 



S 



AtxG 



where Ah represents the area of the black hole and G denotes the Newtonian constant. 

Before proceeding further, a few comments are in order. Since the underlying theory is conformal, 
the only relevant parameter is the ratio T/fi. Thus we define 

Xd = i (~) , (9) 



4tt \T 

which we will consider throughout the paper. 
Now, we would like to convert the metric into Poincare coordinates. By taking 



dv = g(t, z)dt — —. rdz , (10) 

f{v,z) 



where g(t,z) is a hitherto unknown function, the metric in EF coordinates can be rewritten as 



1 / dz 2 
ds 2 = — ( -f(v, z)g(t, z) 2 dt 2 + — r + dx 2 



f(v,z) 



(11) 



The function g(t, z) can now be identified with the redshift factor. In the thin-shell limit we can 



37| and therefore separate the space-time into 



ignore the compression of the shell as indicated in 
the following two regimes: exterior and interior. For the exterior geometry, g(t, z < zq) = 1 where 
Zq denotes the position of the shell. This is because the exterior of the shell is described by an AdS- 
RN geometry. For the interior of the shell, g(t, z > z ) = f(v > v ,z ), where f{v,z) in different 
dimensions are shown below (OQ) and (EJ). We will refer to this as the quasi-AdS (qAdS) background. 
Thus we get: 



ds 2 = < 



4r (-F(z)dt 2 + |gy + dx 2 ) if v > (z < zq) 
\ (-F(z ) 2 dt 2 + dz 2 + dx 2 ) if v < (z > z ) 



(12) 



where 



F(z) = { 



1 - Mz 3 + (1/2)M 4 / 3 QV for d = 3 , 
1 -Mz 4 + (2/3)M 3 / 2 QV for d = A. 



(13) 



Since the falling velocities of the upper and lower surfaces are the same in the thin-shell limit, the 
position of the falling shell is given by 



t{z ) 



o F(z 



(14) 



Given that the shell never coincides with the future horizon exactly in the Poincare coordinates, we 
approximate the thermalization time as r = t(^o)|z =o.99^; which is shown in FigHJand Figf2]with 
different values of chemical potential measured in the unit of temperature for d = 3 and d = 4, 



respectively. When increasing the chemical potential, the thermalization time decreases. We should 
comment here that the specific choice of Zq = 0.99^ does not affect the qualitative feature as far as 
the thermalization time is concerned, any other choice yields a similar behavior. Also, for very large 
value of Xd the thermalization time vanishes asymptotically. 
The qualitative behaviors here are distinct from those found by studying non-local operators in 



the AdS-RN-Vaidya metric [35|. However, we should notice that the definitions of the thermalization 
time in these two approaches are different. For example, if we introduce a spacelike geodesic with 
two ends fixed on the boundary, the thermalization time can be defined as the time when the shell 
grazes the bottom of the geodesic in the bulk. After that, the whole geodesic will reside in the 
AdS-RN spacetime, which matches the result derived in the thermalized medium. Therefore, in 
this situation, the thermalization time of non-local operators are recorded by the position of the 
falling shell. Nevertheless, when the separation on the boundary of the spacelike geodesic is too 
large such that the bottom of it penetrates the future horizon, the geodesic can never fully reside 
in the AdS-RN spacetime even when the shell asymptotically reaches the future horizon. In other 
words, the maximum thermalization time in Poincare coordinates should be the time when the shell 
almost reaches the future horizon. For the non-local operators with much larger lengths can never 
thermalize in Poincare coordinates or the definition of the thermalization time in such a case should 
be modified. Thus, thermalization times extracted from two approaches can only be compared when 
the length scale of non-local operators introduced in EF coordinates is about the size of the future 
horizon. 
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FIG. 1: The thermalization time r with different values of FIG. 2: The thermalization time r with different values of 
chemical potential in d — 3. chemical potential in d = 4. 
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III. THE JET QUENCHING OF VIRTUAL GLUONS 

Let us also study thermalization of light probes traversing the non-equilibrium medium. In this 
section, we will follow the approach in 44J to compute the maximum stopping distance led by the 
falling massless particle characterizing the tip of a double string with two ends fixed at an infrared 
scale in the 5-dimensional AdS-RN-Vaidya spacetime. In the gauge theory dual, this double string 
corresponds to a virtual gluon, in which its energy is encoded by the position of the tip of the string. 
The computation in the 4-dimensional spacetime can be carried out in the same manner. A similar 
computation investigating the stopping distance of a massless particle propagating in the bulk as the 
holographic dual of a R-charged current diffusing on the boundary and the falling string scenario in 



the AdS-Vaidya spacetime is shown in 



37] . 



A. Stopping distance in the AdS-RN spacetime 

Before proceeding with the non-equilibrium medium, we may begin by investigating the thermal 
medium with a non-zero chemical potential, which is described by the AdS-RN geometry. When 
the particle moves in the spacetime metric which only depends on the z coordinate, the energy anc 



spatial momentum of the particle should be conserved. The null geodesic can thus be written as 

dx % gVqj 



dz {-qkqig kl ) 1/2 



47] 



(15) 



where gij is the spacetime metric and qi = (—oj, \q\, 0, 0) for i — 0, 1, 2, 3 is defined as the conserved 
4-momentum of the particle. Here we define the momentum of the particle by using q^ = g^ v dx v /d\, 
A being the affine parameter, in which case we can write down q z in terms of q^ and z. In this note, 
we will use Roman indices and Greek indices to represent 4-dimensional and 5-dimensional spacetime 



9 

directions, respectively, unless otherwise mentioned. In the AdS-RN spacetime, ( !T5|) leads to 

dx 1 1 



(16) 



dt 1 



The first equation in ([16]) allows us to study the stopping distance of the light probe in the thermalized 
medium with a non-zero chemical potential. The stopping distance is given by 

*i = r - — , (17) 

which only depends on the ratio \q\/oJ and the initial position of the tip of the string zj. The ratio 
\q\/oo roughly represents the ratio to the spatial momentum and energy of the virtual gluon and Zj can 
be associated with its initial energy. More comprehensive discussions over the physical interpretation 



of these two quantities can be found in 37j. Even though determining the ratio of |<f| to u; entails the 



initial string profile close to the tip, which could depend on the geometry of the thermalized medium 



38 



the profile of a moving 



we may choose the straight string as the simplest setup As shown in 
string obtained from extremizing the Nambu-Goto action in the AdS geometry is a straight string. 
Although the physical solution in the AdS-Schwarzschild geometry corresponds to a trailing profile, 
where the profile in the presence of a chemical potential could be more complicated, we will make 
the same straight-string setup in the thermalized case for comparison. This assures that the gluons 
traveling in the thermalized medium with different values of chemical potential carry the same initial 
energy and momentum and as well the same ratio to |g| and u. The ratio of the stopping distance 
with a non-zero chemical potential to that with zero chemical potential at the same temperature is 
shown in Figj3] and FigJD for both d = 3 and d — 4. It turns out that the stopping distance of the 
light probe decreases when the chemical potential increases, which has the same qualitative feature 
compared to the thermalization of the medium. 
Physically, by increasing the chemical potential, we actually increase the number of states. From 



10 

(JSJ), the entropy density of the plasma in d + 1 dimension is given by AGs = z~^ d ~ l \ As illustrated 
in Figj5] and Figj6j the entropy density increases when the chemical potential is increased. The 
medium thus becomes denser, which results in the enhanced scattering for the light probe and a 
smaller stopping distance. On the other hand, the medium also thermalizes faster because of the 
same effect. Although the holographic correspondence of the falling shell in the gauge theory side is 
unknown, the falling shell could be characterized by the collective motion of massless particles with 
large virtuality in the gravity dual since a null shell is homogeneous along spatial directions and it 
falls along a null geodesic. Each particle as a component of the shell may as well be influenced by 
the backreaction to the spacetime metric caused by the falling shell; hence it should qualitatively 
behave in the same manner as the light probe traversing the medium. Since the component particles 
of the shell carry large virtuality, the enhanced scattering led by increasing chemical potential would 
be more pronounced, which reduces the thermalization time of the medium. By comparing Figj5] 
with FigJHl the entropy density for d = 4 increases more rapidly than that for d = 3, which also 
manifests the steeper drop of the thermalization time of both the probe and the medium for d = 4 
when increasing the chemical potential by comparing FigJT] with Figj2] and Figj3] with FigJH 

To avoid a naked singularity in the AdS-RN spacetime, there exists a maximum value of Q such 
that the position of the horizon is given by F(z) = 0. By taking M = 1, the maximum values are 
Qmax = (27/32) 1//6 and Qmax = 3~ 1//4 for d = 3 and d = 4, respectively. Given that the derivative of 
F(z) vanishes at the horizon, the temperature of the medium then vanishes. The corresponding zero 
temperature black hole is known as the extremal black hole. However, this does not correspond to 
the vacuum since it carries a non-zero entropy given by a non- vanishing area of the extremal black 
hole horizon. In this situation, we see that the dimensionless parameter X4(3) diverges. Also, the 
stopping distances with a non- zero chemical potential in the unit of temperature, xl(T,fi) shown in 
Figl3]and FigJH become zero; we hence lose the ability to make a comparison between the observables 
in the media with and without a chemical potential in this particular situation. 
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FIG. 3: The ratio to the stopping distances with and without FIG. 4: The ratio to the stopping distances with and without 
chemical potential in the unit of temperature for d = 3, chemical potential in the unit of temperature for d = 4, 
where x s = x\T. Here we set M = 1, zi = 0, and |g| = where x\ = x\T. Here we set M = 1, Zi = 0, and \q\ = 
0.99u. 0.99oj. 
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FIG. 5: The entropy density with different values of the FIG. 6: The entropy density with different values of the 
chemical potential for d — 3. Here we set M = 1. chemical potential for d = 4. Here we set M — 1. 

B. Initial conditions of the falling string and the matching condition 

Now we discuss the physics when the system is not in equilibrium. In the qAdS spacetime, energy 
conservation no longer holds since the redshift factor encoded in the spacetime metric depends on 
the position of the falling shell. However, by redefining dt = F(zo)dt, the spacetime metric is 
independent of the new time coordinate t. By taking the conserved 4-momentum inside the shell 
cji = (—a), \q\, 0, 0), the null geodesic in the quasi-AdS spacetime can be written as 



dx 1 
dz 

dt 



1 



l<?T 



- 1 



1/2 ' 
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Now u and u are two conserved energies denned in the AdS-RN spacetime and the quasi-AdS 
spacetime, respectively. 

When the string fully resides in the quasi-AdS spacetime, the string profile obtained from extrem- 
izing the Nambu-Goto action by choo sing the redshift time coordinate t and z as the worldsheet 



coordinates should be a straight string 37| : x l = vt. The momentum density and momentum of the 
corresponding string take the form: 



f°° i 

p, = 2 dz*l = - (-1, v, 0, 0, 0) , (19) 

J zi ira'zjy/l — v z 

where x M = (dx^/dt) and zj represents the initial position of the tip of the string. It is argued in 



371 ] that the 4-momentum of the massless particle should be proportional to the 4-momentum of 
the string since the falling string here behaves approximately as a null string. We thus have the 
following relation, 



|i = ?i = -v. (20) 

<?o Po 



As pointed out in [37|, the part of a falling string outside the shell could be distorted and bears an 
unknown profile, while the part inside the shell should remain straight. Nevertheless, the detailed 
structure of the string profile should not affect the maximum stopping distance. 

Next, when the tip of the string penetrates the shell and falls in the AdS-RN geometry, we have 
to find the matching condition to connect the u and u defined in two spacetimes at the collision 



point z c . As shown in 



37| , such a condition can be derived by requiring the momentum conservation 



along the spacetime direction tangent to the shell, which results in 



Co' 



W 2(1 



[2F(z c )(l - VT^P) + 5 2 (1 - F{z c )% , (21) 



where 5 = \q\/oo. The ratios are illustrated in FigJ7]and FigJHlfor both d = 3 and d = 4, respectively. 
The same matching condition can also be derived in a different setup by using the continuity of the 



wave function of the massless particle in the WKB approximation 37 ]. 
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FIG. 7: The red, blue, and green solid curves represent the FIG. 8: The red, blue, and green solid curves represent the 
energy ratios for Q = with |<fj = 0.99, 0.9, and 0.5u;, re- energy ratios for Q — with \q\ = 0.99, 0.9, and 0.5tl>, re- 
spectively. The case for Q = 0.9 are showed by the dashed spectively. The case for Q = 0.7 are showed by the dashed 
curves. Here we set M = 1 and d = 3 curves. Here we set M = 1 and d = 4 



C. Stopping distance in the AdS-RN-Vaidya spacetime 



Now, we may follow the approach in 37J to compute the stopping distance in the thin-shell limit 
in Poincare coordinates. We will also eject the shell and massless particle at the same time. The 
computation involves tracking the shell and the massless particle simultaneously and finding the first 
collision point z c . For simplicity, here we only show the general expression for the stopping distance, 

! _ f Zc dz f Zb dz 

X s — I , „ 1/2 + / f v 1/2 ' 



'w- 1 ) Jzc [w- F V< 



where the second collision point z^ is approximately equal to Zh since the shell falls with the speed of 



37] . the stopping distance of the massless 



light. Similar to the study in the AdS-Vaidya spacetime 
particle falling from the boundary is not affected by the gravitational collapse in the AdS-RN-Vaidya 
geometry In our setup, when we set the initial position of the tip of the string on the boundary, 
the corresponding virtual gluon on the gauge theory side carries infinite energy as shown in f )19p . As 



indicated in [37J, when the hard probe has the energy much larger than any other scale such as the 
thermalization temperature or the chemical potential of the system, the thermalization of the probe 
would be insensitive to the thermalization of the medium. However, for the soft probe carrying 
energy comparable to other scales of the system, we may envision an influence of the thermalization 
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371 ] for a zero chemical potential. 



process on the jet quenching phenomenon, which is also explored in 
This scenario is shown in Fig JS] and Figfiuj where we eject the massless particle below the boundary, 
which corresponds to the soft gluon with a finite energy. As a result of thermalization, stopping 
distance of light probes increases. 

Finally, we illustrate this behaviour for different values of chemical potential in AdS-RN and 
AdS-RN- Vaidya spacetimes in FigJTT] and Fig{12j The results can be reproduced by solving geodesic 
equations directly in EF coordinates, which are presented in Appendix A; the latter approach can be 
further applied beyond the thin-shell limit. As shown in FigJTTJand FigJT2l the stopping distances 
scaled by the thermalization temperature in both the AdS-RN- Vaidya and AdS-RN geometries 
decrease when the values of \d are increased and will drop to zero when the values of \d reach 
infinity. 

In general, we find that the probe gluon travels further in the non-equilibrium plasma with a 
non-zero chemical potential compared to the thermal background. Increasing the magnitude of the 
chemical potential decreases the stopping distances in both equilibrium and non-equilibrium plasmas. 
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FIG. 9: The red and blue curves represent the trajectories 
of the massless particles moving in AdS-RN and AdS-RN- 
Vaidya spacetimes for d = 3 and X3 = 4.47, respectively. 
The red and blue dashed lines denote the first collision point 
and the position of the future horizon. Here we take zi — 0.4, 
M — 1, and \q\/& = 0.99 as the initial conditions in both 
spacetimes. 



FIG. 10: The red and blue curves represent the trajectories 
of the massless particles moving in AdS-RN and AdS-RN- 
Vaidya spacetimes for d = 4 and \± = 1-1, respectively. The 
red and blue dashed lines denote the first collision point and 
the position of the future horizon. Here we take zi = 0.4, 
M — 1, and \q\/oJ = 0.99 as the initial conditions in both 
spacetimes. 
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FIG. 11: The red and blue points represent the stopping FIG. 12: The red and blue points represent the stopping 
distances with different values of chemical potential in AdS- distances with different values of chemical potential in AdS- 
RN and AdS-RN-Vaidya for d = 3, respectively. Here we set RN and AdS-RN-Vaidya for d = 4, respectively. Here we set 
M = 1, Zi = 0.4, and |g|/tl> = 0.99 as the initial conditions M = 1, Zi = 0.4, and \q\/u> = 0.99 as the initial conditions 
in both spacetimes. in both spacetimes. 



IV. SUMMARY AND DISCUSSIONS 



In this paper we have investigated thermalization of a non-equilibrium plasma with a non-zero 
chemical potential by tracking a thin shell falling in the AdS-RN-Vaidya spacetime. We found the 
thermalization time of the medium decreases when the chemical potential is increased. We have also 
studied the jet quenching of a virtual gluon traversing such a medium by computing the stopping 
distance of a falling string, in which the tip of the string falls along a null geodesic. In both the 
thermalized or thermalizing medium with a non-zero chemical potential, the stopping distance of 
the probe gluon decreases when the chemical potential is increased. On the other hand, for a soft 
gluon with finite energy comparable to the thermalization temperature of the medium, its stopping 
distance in the thermalizing medium is larger than that in the thermalized case. 

In section II, we briefly discussed the difference between the thermalization time obtained from 
our approach and that derived in 3s| ; we will further elaborate on this here. When the medium 
carries no chemical potential, the position of the horizon is about the inverse of the temperature. 



As indicated in [37|, the thermalization times obtained from two approaches in the AdS-Vaidya 
geometry approximately match when the length scale of the nonlocal operators is about the inverse 
of the temperature. However, for the medium with a non-zero chemical potential, the temperature 
does not linearly depend on z^ 1 . To compare the thermalization times obtained from two approaches, 
we have to investigate the thermlization time for a nonlocal operator with the length scale about 
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the size of the horizon in the AdS-RN-Vaidya spacetime. As shown in Fig JT3l t hermalization time 

35[ do exhibit similar 



obtained by analyzing non-local observables with a length-scale l s = \.l\Zh 
qualitative feature as the one we have encountered here. Note that the number l s j Zh bears no possible 
physical significance other than being an order one number where we have found a visibly pleasant 
matching. As far as the thermalization time of the medium is concerned, we conclude that our 
approach by tracking the falling shell close to the horizon seems consistent with that by probing the 
thermalizing medium with nonlocal operators, when the length scale of the operators approximately 
equals the size of the horizon. In addition, the thermalization times from both approaches starts to 
drop more rapidly when [i > T. 

For the light probe traversing the medium with a non-zero chemical potential, the decrease of the 
stopping distance when increasing the chemical potential is expected due to the enhanced scattering 



with the increasing density of the medium. In 



55 



561 ]. it is found that the jet quenching parameter 



and the drag force of a trailing string in the charged SYM plasma both increase when the chemical 
potential is increased, which is again consistent with our general observations here. 




0.25 
0.20 



0.2 0.4 0.6 0.8 1.0 



FIG. 13: The blue and red points represent the thermalization times scaled by the temperature obtained from our approach 
and that from analyzing non-local observables, respectively. Here we take M = 1 in the both cases. 
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VI. APPENDIX 

A. Finding the stopping distance in Eddington-Finkelstein coordinates 

It is technically difficult to track the null geodesic of the massless particle beyond the thin shell limit 
in Poincare coordinates, since the trajectory of the particle is continuously deflected by interactions 
with the shell. However, irrespective of the details of the interaction, it is easier to compute the 
stopping distance of the probe in EF coordinates. A similar computation in the AdS-Vaidya geometry 



has been done in [37j]. By using the metric in ([I]) or (j2J) and setting the momentum component of 
the massless particle non-zero only along v, z, and x 1 directions, we can write down the geodesic 
equations in terms of the affine parameter A, 

d 2 v „„, / dv \ 2 „„, ( dx 1 ^ 



d\ 2 vv \dXJ u \dX 



^l l (^)(^)=Q. (23) 



d\ 2 vv \dX J vz \dX I \dX J \ dX J zz \dX 

dX 2 zl \dX J \dX 

For d — 3, the relevant Chiristoffel symbols are given by 



K v = - (-1 + \m{v)z* + \m{vfz 6 - jQ 2 m(v) 7 ' 3 z 7 + l -z*Q 4 m(vf" 



1 f 1 ^ i^/.Ai/a A dm{v) 



ZZ 1 

Z 



T z vz = -i (l + m(v)z 3 - l -Q 2 m{vY"z^ , V 
rfi = ; (l - ™(^ 3 " ^Q 2 m(^) 4/3 ^ 4 ) , I* = -i . 
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For d — 4, 



Kv = ~ (l + m(v)z* - -Q 2 m{vfl 2 z") , V u = -- , 
z \ 6 J z 

T z vv = -(-1 + m(v) 2 z 8 + -Q 2 m(vf 2 z 6 - 2Q 2 m(vf 2 z w + -Q 4 m(v) 3 z 
z \ 3 9 



2z ' dv 



(25) 



\Q 2 m{vfl 2 z^ , 

r?i = - f i - m( v y + 2 -Q 2 m(v f 2 z«) , i* = -1 . 

2; \ 6 J z 

In addition to these, the mass function is defined as 

m W = f (l + tanhg-i)) , (26) 

where the shift in the hyperbolic-tangent function is to fit the setup illustrated in FigfHl In the 
thin-shell limit, the mass function defined here will degenerate to the expression in (j4]). The mo- 
mentum along the x 1 direction q± = \q\ = z ~ 2 ^ from the definition = g^u^ is conserved in EF 
coordinates. Thus, we may rewrite the geodesic equations in terms of the derivative with respect to 
x 1 , 

dx^ dx^ dx 2 1 — *i 1 ^1 
~d\ = dx^l\ =Z |g> ' 

= \q\ 2 z 2 (x"»z 2 + 2zz>x>n , (27) 

where the prime denotes the derivative with respect to x . By inserting ( 1271) into the geodesic 
equations, we find the last equation in ( 1231) is automatically satisfied due to the conservation of 
momentum along x 1 . We then have to introduce proper initial conditions to solve other two geodesic 
equations in (1231) . Even though the energy of the particle is not conserved in the AdS-RN-Vaidya 
spacetime, initially the gravitational effect led by the falling shell is negligible. Therefore, we may 
write down the initial conditions as a massless particle falling in the pure AdS metric. For a particle 
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ejected from zj, we have 



v Ui=o 



^ a;l=0 



\x 1 =Q 



— Z\~i 



J-!=0 



s 1 =0 



\q\ 



{21 



where uq and q represent the initial energy and momentum of the massless particle in Poincare 
coordinates. The stopping distance obtained in EF coordinates is shown in Fig|T5]for d = 4 in the 
thin-shell limit, where we also make a comparison with the result illustrated in Fig{12]with the same 
initial conditions. The results match and cannot be distinguished from the plot. This holds for d = 3 
as well. We can also evaluate the stopping distance with the thick shell in EF coordinates, it turns 

out that the deviation from the thin shell is negligible. 

z 



shell 




s - 

f- — 



0.5 
0.4 
0.3 
0.2 
0.1 



0.2 0.4 0.6 0.8 



1.0 



X4 



FIG. 14: The scenario of a massless particle electing from „ T „ , _ . , , , , ,. . 

, , , , , , , „ . „ , . 11G. 15: 1 he sreen points represent the stopping distances m 

the boundary as the shell starts to rail, where vn denotes . ,„ _. AT , . . , .. . „„ ,. , , . , . , 

, ... . . ... , . .. . , , , Adb-KN-vaidya spacctime m lil coordinates, which match 

the thickness ot the shell and the solid red curves and the , . ,. , , .. . , 

.... , , , r , , . . , . those derived m Poincare coordinates as shown by the blue 

dashed red curve represent the surfaces and the center ol the . , . . ,, ,, 

, , ■ , r t iiii- ,i points. Here the initial conditions are the same as those m 

shell, respectively. Here the blue dashed line represents the zi. — „ . , . 

. j_"i j_ i n FiglT2land we take vn = 0.0001. 

masselss particle ejected irom zi. 
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B. The dyonic black hole 

Here we will summarize some properties of the dyonic black hole in d = 3. The action we consider 
is the following 



(29) 



which results in the following eom 

±(R-2A)g^ = g a "F m F av -± 
d p [y/=jF"]=0, 0, [V=£(*iTj = , (31) 



- - (R - 2A) 9iiV = g ap F pfl F au - -g^F 2 , (30) 



where (*F) is the Hodge dual of F. In (3 + l)-bulk dimensions, the Hodge dual of the bulk electro- 
magnetic field is another two form. Thus it is possible for a black hole to possess both electric and 
magnetic charges. Such a dyonic black hole is given by 

ds 2 = ^(-fdt 2 + ^ + dx 2 ^J , A = — , (32) 

f( z ) = l-Mz s + ^(Ql + Q 2 m )z 4 , (33) 
F = Q e dz A dt + Q m dx A dy , (34) 

where Q e and Q m denote the electric and the magnetic charges respectively. 
In the extremal case, the function / can be written as 

f(z) = l-J-X + 3(-X , (35) 



\ Z H / \ZH 

where Zh is the location of the event-horizon. It is easy to check that the temperature is given by 



4vrT = -# 

dz 



= . (36) 
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Therefore for the extremal case, we can write 



Z ^ ~ M ' ^ + Qm — j2 „4 • (^) 



We can now find the corresponding Vaidya-type background sourced by appropriate matter field 

S = S + K^cxt • (38) 

This will give the following equations of motion 

Rp, -\i R - 2A ) 9,u - 9 ap F Pil F au + ± 9lu ,F 2 = (16ttG n k) , (39) 

d p [yf^F^] = (8ttG n k) r e , (40) 

d p [V=g{*F)"] = (8ttG n k) r m . (41) 

Here J e denotes the electric current and J rn denotes the magnetic current. The dyonic-Vaidya 
background takes the following form 

ds 2 = ^ (-fdv 2 - 2dvdz + dx 2 ) , f(z, v) = l- m{v)z 3 + ^ (q e (v) 2 + q m {vf) z A , (42) 

with the following vector fields 

Fzv = q e (v) , F xy = q m (v) . (43) 
The above background is sourced by the following stress-energy tensor 

«f = 4 - * (?«% + J£S\ S m S„ . (44) 



The electric and magnetic sources are given by 



K J£ = <^± 5 >™ kJ m = J^n^v (45) 

dv dv 
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The extremal limit is now obtained by considering 




4 



Qe{v) 2 + q m {v) 2 



6 1 



(46) 




Once m(v) is chosen, we can pick the electric and magnetic charge functions according to the above 
formula. An obvious such choice is given by 
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